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ABSTRACT 

The  rheological  properties  of  two  sulfur  -  water  suspensions 
with  approximately  the  same  particle  size  analysis  have  been  investigated 
using  a  modified  Fann  concentric  cylinder  viscometer.  Both  suspensions 
gave  flow  curves  similar  to  those  expected  for  pseudo plastics.  The 
experimental  data  are  presented  as  logorithmic  plots  of  shear  stress 
versus  shear  rate. 

Pressure  drop  -  flow  rate  data  have  been  obtained  over  a 
range  of  solids  concentrations  for  suspensions  flowing  through  40.17 
foot  and  45*35  foot  lengths  of  1-7/8  I.D.  smooth  aluminum  tubing.  The 
pressure  drop  data  were  obtained  using  an  electrically  operated  differ¬ 
ential  pressure  cell  associated  with  a  potentiometer. 

A  comparison  has  been  made  between  the  experimental  results 
and  the  predicted  values  assuming  the  theory  applicable  both  for 
Bingham-plastic  behaviour  and  for  "power-law”  pseudoplastic  behaviour. 

As  a  result  of  these  comparisons  it  is  concluded  that  one 
of  the  suspensions  can  be  treated  as  a  Bingham-plastic  in  the  range 
of  the  pipeline  data.  The  properties  of  the  other  suspension  are 
more  complex,  it  behaves  as  a  pseudoplastic  that  does  not  obey  the 


power  law. 
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INTRODUCTION 

During  the  last  few  years  the  natural  gas  industry  in  the 
Province  of  Alberta  has  developed  rapidly.  A  by-product  from  this 
industry  is  sulphur  and  this  is  now  being  produced  in  large  quantities. 
In  order  to  avoid  an  accumulation  of  large  quantities  of  this  sulphur 
a  method  of  transportation  will  have  to  be  found.  A  convenient  method 
of  transporting  this  material  would  be  by  pipeline  in  the  form  of 
either  a  molten  liquid  or  a  solid-liquid  suspension. 

The  determination  of  commercially  useful  data  was  not  the 
aim  of  this  investigation.  In  the  initial  stages  it  was  planned  to 
combine  the  obtaining  of  rheological  and  pressure  drop  data  with  an 
attempt  to  correlate  the  effects  of  particle  size  on  the  pressure 
drop,  however,  due  to  inadequate  equipment,  this  was  not  accomplished. 
Instead  a  comparison  has  been  made  between  the  experimental  results 
and  values  predicted  by  theoretical  equations. 
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LITERATURE  SURVEY  AM)  THEORY 
I  Fluid  Classifications, 

Fluids  are  classified  as  Newtonian  or  non-Newtonian 
according  to  their  behavior  at  constant  temperature  under  imposed 
shearing  forces  in  laminar  flow.  To  illustrate  this  behavior 
suppose  a  fluid  is  confined  between  two  plates  of  area  A,  and 
separated  by  a  differential  distance  dr.  If  the  lower  plate  is 
fixed,  a  small  force  F  applied  to  the  upper  plate  will  give  it  a 
differential  velocity  dV  in  the  direction  of  the  force.  A  uniform 
velocity  gradient  of  magnitude  dV/dr  is  produced  in  the  fluid  since 
the  shearing  force  is  uniform  across  the  distance  dr. 

The  velocity  gradient  dV/dr  is  referred  to  as  the  "rate  of 
shear",  similarly,  the  shearing  force  per  -unit  area,  F/A,  is  known 
as  the  "unit  shearing  stress"  and  is  denoted  by  T. 
a)  Newtonian  Fluids: 

By  definition  these  fluids  follow  the  equation 


dr 

Thus,  for  laminar  flow,  there  is  a  direct  proportionality  between 
shear  stress  and  shear  rate.  The  constant  of  proportionality,  p,  is 
known  as  the  viscosity  of  the  fluid. 

This  equation  is  illustrated  graphically  in  Figures  1  and  2, 
where  shear  stress  is  plotted  against  shear  rate  on  arithmetic  and 
logarithmic  co-ordinates  respectively.  On  arithmetic  co-ordinates 


(1) 
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the  equation  is  represented  by  a  straight  line  of  slope  p  passing 
through  the  origin.  On  logarithmic  co-ordinates  the  curves  for  all 
Newtonian  fluids  have  a  slope  of  unity  while  the  viscosity  is  given 
by  the  intercept  at  a  shear  rate  of  unity, 
b)  Non-Newtonian  Fluids: 

Non-Newtonian  fluids  do  not  obey  Equation  1  and  so  do  not 
possess  the  property  of  viscosity.  The  measure  of  the  consistency  of 
these  fluids  is  known  as  either  the  "apparent  viscosity"  (7)  or  the 
"coefficient  of  rigidity"  (17). 

There  are  five  general  types  of  non-Newtonian  fluid 

1)  Bingham  Plastic 

2)  Pseudoplastic 

3 )  Dilatant 

4 )  Thixotropic 

5)  Rheopectie 

the  latter  two  types  being  more  complex  than  the  first  three  in  that 
the  relationship  between  shear  stress  and  shear  rate  depends  upon  the 
duration  of  shear . 
l)  Bingham  Plastic  Fluids: 

As  shown  in  Figure  1  this  is  the  non-Newtonian  which  differs 
least  from  a  Newtonian  fluid.  The  relation  between  shear  stress  and 
shear  rate  is  still  linear  but  it  does  not  pass  through  the  origin. 

Thus  a  finite  shear  stress  Ty,  known  as  the  yield  stress,  is  necessary 
to  cause  motion.  The  equation  to  the  flow  curve  becomes 

T  -  Ty  =Y)  JX  . (2) 

^  dr 
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where  yj^  ,  the  "coefficient  of  rigidity",  or  "apparent  viscosity", 
represents  the  slope  of  the  flow  curve.  On  logarithmic  co-ordinates 
the  flow  curve  becomes  asymptotic  to  Ty  at  low  shear  rates,  and 
approaches  a  slope  of  unity  at  very  high  shear  rates. 

2)  Psuedoplastic  Fluids: 

On  arithmetic  co-ordinates  a  pseudoplastic  fluid  displays 
a  flow  curve  which  is  concave  to  the  rate  of  shear  axis  and  passes 
through  the  origin  (Figure  l).  When  represented  on  logarithmic 
co-ordinates  (Figure  2)  these  fluids  exhibit  flow  curves  having  slopes 
between  zero  and  unity,  and  these  are  generally  relatively  straight 
over  wide  ranges  of  shear  rates  (l6).  Certain  pseudoplastic  fluids 
exhibit  flow  curves  which  can  be  represented  by  the  power  law  function 


T  =  K  (-  ^)n 
dr 


These  materials  are  known  as  "power-law"  pseudoplastics. 

3)  Dilatant  Fluids: 

Dilatant  fluids  display  rheological  behavior  opposite  to 
that  of  pseudoplastics  (Figures  1  and  2)  in  that  the  flow  curves  are 
convex  to  the  rate  of  shear  axis. 
k)  Thixotropic  Fluids: 

When  thixotropic  fluids  are  sheared  at  a  given  shear  rate 
and  temperature  the  shear  stress  decreases  with  the  duration  cf  shear. 
5)  Rheopectic  Fluids: 

These  fluids  behave  in  the  converse  manner  to  thixotropic 
fluids  in  that  the  shear  stress  increases  with  the  duration  of  shear 


(3) 


at  constant  shear  rate. 
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Examples  of  each  of  these  fluids  and  explanations  of  the 
causes  of  their  rheological  behaviour  are  available  (l,  10,  l4,  17, 
25,  26). 

II  Mechanics  of  Flow  in  Circular  Pipes. 

Relations  used  for  the  flow  of  fluids  are  developed  from 
material  balances,  force  balances  and  energy  balances  (24). 

Applying  a  force  balance  to  steady  uni-dimensional  flow  of 
any  fluid  in  a  pipe  yields 


V  dV  4  T-w  v  dL  g  da 

-  v  dp  =  -  +  - “ -  +  — 


a  gr 


D 


(4) 


Jc  c 

A  non-dimensional  shear  stress,  f ,  known  as  the  friction  factor,  is 
defined  by  the  relationship 

2  gc  ^w 


f  = 


P 


(5) 


By  dimensional  analysis  this  friction  factor  may  be  shown  to  be  a  function 
of  the  relative  roughness,  6/D,  of  the  pipe  surface,  and  the  Reynolds 
number,  DVp/(a. 

Eliminating  T  from  Equations  3  and  4  gives 

V  dV  2f  V2  dL  g  da 

-dp-  +  +  .o... 

a  v  gc  gc  D  v  gc  v 

these  three  terms  represent  the  loss  in  pressure  due  to  acceleration, 
friction  and  change  in  elevation  respectively,  so 


dpf 


2  f  V2  dL 

gc  D  v 


(7) 


this  is  known  as  the  Fanning  equation 
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Laminar  Flow  in  Circular  Pipes* 

Consider  the  steady  laminar  flow  of  a  Newtonian  fluid  in 
a  horizontal  pipe  of  uniform  diameter.  The  velocity  at  any  given 
distance  r  from  the  axis  of  the  pipe  can  he  obtained  by  a  force 
balance,  and  is  given  by 


(8) 


However  the  volumetric  rate  of  flow  through  an  annulus  of  radii  r 
and  r  +  dr  is 


(9) 


dQ  =  2jrr  dr  ur 


eliminating  ur  from  Equations  8  and  9  gives 


D 


integrating  this  equation  between  the  limits  r  =  ^  and  r  =  0  gives 


Q 


it  D4  dpf 


12511  dL 


Thus  the  average  velocity, 


V 


Q 


V 


D 2  dPf 
32  \i  dL 


or  -  dpf 


32  it  V  dL 


(10) 


This  is  the  well  known  Boiseuille  equation. 
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Now  equating  the  pressure  drop  term  of  Equation  10  with  Equation  7  gives 


f  = 


16 

DV  p 


16 

Re 


(ii) 


i.e.  the  friction  factor  is  proportional  to  the  Reynolds  number  only 
in  laminar  flow. 

Turbulent  Flow  in  Circular  Pipes 

Turbulent  flow  is  a  much  more  complex  condition  than  laminar 
flow  and  is  not  as  easily  analyzed.  However  there  are  several  empirical 
equations  for  the  velocity  profile  which  appear  to  fit  the  data  over  a 
large  percent  of  the  pipe  diameter. 

The  simplest  of  these  equations  is  (24) 

-  u 


U 


max 


-ai£ 


V, 


•  (12) 


*  \  W' 

where  A  has  a  value  of  7* *85.  It  is  found  that  this  equation  does  not 
apply  near  the  pipe  wall. 

The  other  common  expression  for  the  velocity  profile  is  (24) 


11 


max  u  r 

—  =  2.  5  In  — 


V 


(13) 


Lw 


this  equation  also  does  not  apply  very  close  to  the  pipe  wall. 

The  Flow  of  non-Newtonian  Fluids  in  Pipes 

The  literature  relating  to  the  flow  of  non-Newtonian  fluids 
in  circular  pipes  was  reviewed  by  Winning  (29)  in  his  M.Sc»  thesis  in 
1948.  Ten  years  later,  in  1958,  Lindley  (l4)  again  reviewed  the 
literature  and  found  there  had  been  a  considerable  amount  of  new  material 
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published.  Most  notable  of  this  was  a  survey  by  Metzner  (17)  on  the 
whole  field  of  non-Newtonian  technology,  giving  a  detailed  report  on 
the  theory  to  that  date.  Since  then,  further  work  has  been  reported 
especially  related  to  tne  theory  of  the  turbulent  flow  of  non-Newtonian 
fluids  in  pipes, 
i)  Bingham  Plastic  Fluids: 

Due  to  the  complex  nature  of  turbulent  flow  an  empirical 
approach  is  necessary  in  the  development  of  a  relationship  which  can 
be  applied  to  turbulent  as  well  as  laminar  flow.  A  convenient  empirical 
approach  is  that  of  dimensional  analysis.  Thus,  for  a  Newtonian  fluid 

AP  =  O'  (D,  L,  V,  6,  p,  n,  gc) . (14) 

Similarly  for  a  Bingham  plastic 

AP  =  9"  (D,  L,  V,  6,  p,v^,  gc,  Ty)  . (15) 

Govier  and  Winning  (6)  applied  this  technique  to  the  flow  of  Bingham 
plastics  in  a  horizontal  circular  pipe.  They  rearranged  the  variables 
to  form  four  non-dimensional  groups  as  the  equation 


gc  DAP 
2  V2  L  p 


9" 


D  V  p 

-T" 


which  can  be  rewritten  as 

f  =  e''(Rffi  h  y) 


5  D  Ty  gc 
D  . 


gc  D  AP 

where  the  friction  factor,  f  =  - ^ - 

2  V2  L  p 


(16) 


(16a) 

(17) 


a  modified  Reynolds  number,  P^ 

=  D  Ty  Sc 

VYl 


D  V  p 

n. 


(18) 


a  yield  number,  Y 
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A  different  approach  was  followed  by  Buckingham  (3)*  He 
integrated  the  rheological  equation 


(2) 


for  the  isothermal  flow  of  Bingham  plastics  in  circular  pipes.  The 
resulting  equation  is 


x  41 

8V  _  ^w  ^c 

1 

4 

!z  +  1 

D  ^ 

-L 

"  3 

Tw  3 

lv  _ 

Equation  20  is  usually  modified  in  the  following  manner.  For  the 
shearing  stress  at  the  wall,  the  relationship  is 


(20) 


T 


w 


D  AP 

4l 


(21) 


Substitution  of  Equation  21  in  Equation  20  gives 


8v  _  fc  D  AP 


(22) 


Multiplication  of  both  sides  of  Equation  22  by  /SV^p  and  substitution 
of  the  dimensionless  ratios  defined  by  Equations  17,  18  and  19  yield  the 
Buckingham  equation  in  the  form 


1  =  f  _Y_  1  T' 

Em  16  "  6%  '  3  f3  Em^ 

Division  of  both  sides  of  Equation  23  by  f,  yields 


(23) 


_1 _  =1  1  _Y_  +  1 

f  Em  16  "  6  f  Rm  3  (f  R„,) 1 


(2k) 


which  is  one  of  the  most  convenient  forms  of  the  Buckingham  equation 
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Equation  24  can  be  represented  graphically  in  two  ways, 
either  by  a  logarithmic  plot  of  the  product  f  Rm  versus  Y  or  by  a 
logarithmic  plot  of  f  versus  Rm  with  Y  as  a  parameter.  The  former 
gives  a  single  line  curve  which  is  asymptotic  to  f  Rjh  =  l6  at  low 
values  of  Y,  and  approaches  a  straight  line  for  values  of  Y  greater 
than  1000.  The  latter  yields  a  series  of  parallel  lines  with  slopes 

of  minus  one  and  offset  from  the  Poiseuille  line  ( f  =  —  ,  Y  =  0.) 

\  Re  / 

Govier  (8)  programmed  Equation  24  for  digital  computer 
solution  for  wide  ranges  of  values  of  Rm  and  Y  and  presented  large 
scale  logarithmic  plots  of  these  data  as  f  versus  Rm,  and  f  Rm  versus  Y. 

A  similar  approach  was  used  by  Hedstrom  (l2).  He  presented 
his  results  as  a  plot  of  f  versus  Rj^  but  used  a  different  non-dimensional 
quantity  Y  Rm,  as  parameter.  This  parameter-dias  been  called  the  Hedstrom 
number.  Perkins  and  Glick  (19)  later  rearranged  Equation  20  in  terms  of 
this  non-dimensional  group  to  obtain 


1  =  f  _  _H__  H4 

%  16  6  R,,,2  3fIV0 


where 


H  “  Rra  Y 


Ty  D2  gc  p 

^2 


(25) 


In  turbulent  flow  Equation  24  no  longer  applies  and  the 
importance  of  Y  is  diminished.  If  a  logarithmic  plot  (f  versus  Rm)  is 
made  of  turbulent  data  it  is  found  to  fall  close  to  the  Von  Karman- 
Nikuradse  line  with  Rj^  =  D  V  p/v^  replacing  D  V  p/p  as  the  Reynolds  number. 
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ii)  General  Discussion  Applicable  to  all  Fluids: 

Most  development  related  to  the  pipeline  flow  of 
pseudoplastic  materials  stem  from  relationships  developed  by 
Rabinowitsch  (20)  and  Mooney  (l8).  The  latter  shoved  that  the 
relationship  between  the  shear  stress  at  the  vail  of  the  pipe,  D  AP  , 

Oy  ^ 

and  the  term  —  were  independent  of  the  diameter  of  the  pipe  in 
laminar  flow. 

Alves,  Boucher  and  Pigford  (l)  made  use  of  this  relationship 
between  and  21  in  correlating  data  on  non-Newtonian  suspensions 

and  solutions.  For  laminar  flow  they  suggested  that  a  comparison  of 
pipeline  and  rotational  viscometer  data  could  be  made  on  the  basis  of 
comparing  the  rate  of  shear  and  the  shear  stress  at  the  pipe  wall  with 
the  rate  of  shear  and  the  shear  stress  at  the  inner  cylinder  wall  of 
the  viscometer.  They  also  recommended  a  design  procedure  for  turbulent 
flow  which  was  based  upon  experimentally  determined  viscosities. 

Several  other  relationships  and  procedures  have  been  proposed 
as  useful  in  the  prediction  of  pressure  drop  data: 

l)  An  empirical  power-law  function 


T 


(3) 


The  uses  of  this  function  are  dealt  with  in  the  next  section 


2)  The  Powell-Eyring  equation 


(26) 
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has  been  used  by  investigators  at  the  University  of  Utah.  Christiansen, 
Ryan  and  Stevens  (4)  have  integrated  this  equation  numerically  and  have 
presented  the  resulting  relationships  as  functions  of  non-dimensional 
parameters.  However  it  is  applicable  only  for  laminar  flow. 

3)  For  turbulent  flow  an  apparent  viscosity  at  infinite 
shear  has  been  used  (28)  to  calculate  a  modified  Reynolds  number  and 
hence  to  obtain  a  friction  factor  from  the  conventional  friction  factor 
versus  Reynolds  number  plot, 
iii)  Power- Law  Fluids: 

The  starting  point  for  the  theory  on  power-law  fluids  is  an 
equation  developed  by  Rabinowitsch  (20)  and  Mooney  (l8)  for  the  shear 
rate  at  the  wall  of  a  tube.  The  equation  is 


(_  dV)  =  3  (  §V\  +  1  DAP  d(8 V/D)  .  .  . 

V  dr/w  E  V  D  /  4  TIT  d(DAP/4L) 

Equation  27  was  rearranged  by  Metzner  and  Reed  (l6)  to  give 


(27) 


(-  dV\  = 

l  dr  4 


3n*  +  1 

4n* 


(?) 


(28) 


where  n 1  = 


/  D  AP\ 
d\ln  ~^L~ ) 

* (in  SZ) 


(29) 


D  AP  8V 

which  is  the  slope  of  a  logarithmic  plot  of  ^ ^  versus 


It 


follcn-r  s  that 


D  AP 

4l 


(30) 


r  ; 


i ..V 


» . 
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If  n *  is  a  constant  then  this  equation  is  rigorously  true;  if  n'  is 
not  a  constant  but  varies  with  Tw  then  Equation  30  represents  the 
tangent  to  the  curve  at  any  chosen  value  of  . 

The  similarity  between  Equations  3  and  30  is  obvious; 
however  Equation  30  has  the  advantage  that  it  does  not  have  to  be 
integrated  in  order  to  relate  pressure  drop  to  flow  rate.  Metzner 
and  Reed  (l6)  have  shown  that,  for  fluids  that  obey  the  power- law 
(Equation  3  ) 

n*  =  n  . (31) 


and  K'  =  K 


3n*  +  l\ 

~T~) 


n' 


(32) 


if  the  shear  stress  -  shear  rate  data  has  been  obtained  from  a 

rotational  viscometer.  These  equations  may  be  used  with  Equation  30 

.  .  .  DAP  8V 

to  draw  a  curve  of  ■  versus  —  . 

The  parameters  n  and  n'  characterize  the  extent  of  non- 

Newtonian  behaviour.  When  they  equal  unity,  the  fluid  is  Newtonian; 

values  of  n  and  n*  less  than  -unity  indicate  pseudoplastic  or  Bingham- 

plastic  behaviour,  and  values  greater  than  unity  indicate  dilatant 

behaviour.  The  use  of  exponents  such  as  n  and  n*  has  also  been  proposed 

by  Weltmann  (27). 

Regrouping  Equation  5  as 


D  AP  / 

"  TT/ 


pJL2 

2gc 
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(33) 
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and.  solving  for  a  generalized  or  power-law  Reynolds  number  using 
Equations  11,  30  and  33  yields 

Dn'  V2"n*  « 


Re 


t 


* 


where  V  =  gc  K*  8 


n*-l 


(34) 

(35) 


The  denominator  ^  can  be  obtained  from  rotational  viscometer  data, 
for  power-law  fluids,  by  using  Equations  3 1>  32  and  35. 

For  Newtonian  fluids,  n*  =1,  K*  =  p/gc  and  Equation  34 

reduces  to 

Re-  -  Re  -  2-2-E 

M- 

1  ^ 

As  the  equation  f  =  —  has  been  used  in  the  derivation 

Re 

of  Equation  34  it  follows  that  a  power-law  fluid,  in  laminar  flow, 

should  give  data  that  falls  on  a  line  given  by  f  =  —  ,  when  plotted 

Re’ 

as  friction  factor  versus  power  law  Reynolds  number.  Metzner  (17)  has 
reviewed  a  large  amount  of  published  data  for  the  laminar  flow  of  power- 

law  solutions  and  suspensions  and  has  obtained  good  agreement  with  the 

.  .  .  16 

equation  f  =  — t  . 

A  convenient  approach  to  the  development  of  a  relationship 
for  "power-law1’  pseudoplastics  is  dimensional  analysis.  Thus 

AP  =  ©»*•  (D.  v.  p.  K.  L.  n‘) . (35A) 

Rearrangement  of  these  variables  yields  three  dimensionless  groups  as 
the  equation 


f  Re’  =  2 

which  should  be  applicable  to  laminar  and  turbulent  flow. 
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Ryan  and  Johnson  (22)  have  suggested  a  local  stability 


parameter 


Z 


rw  P  u  du 


(36) 


^  dy 


as  the  starting  point  for  the  estimation  of  the  pressure  drop  in  the 
transition  from  laminar  to  turbulent  flow.  This  parameter  is  applicable 
to  Newtonian  or  non-Newtonian  flow.  Application  of  Equation  36  requires 
knowledge  of  the  laminar  velocity  profile,  which  can  be  derived  from 
the  rheological  law  for  the  fluid  of  interest. 


For  a  Newtonian  fluid 


and  as  the  critical  value  of  the  Reynolds  number  is  2100  then 


(37) 


2100  =  808 


Ryan  and  Johnson  believe  that  808  is  the  critical  value  of  the  parameter 
irrespective  of  the  rheological  behaviour  of  the  fluid. 


For  a  power- law  fluid 


(38) 


at  r 


n/n  +  1 


2 


(39) 


-  IT  - 


By  definition  the  friction  factor  is  given  by  the  expression 


f  = 


2  ^ 

~~5  2 

rw  P  W 


(5) 


Combining  Equations  5  and  38  and  setting  (Zmax)c  =  808  gives 


f  =  A  (a) 

fc  404 


•  o 


(40) 


Experimental  data  may  be  used  to  check  the  accuracy  of  this  equation 
by  plotting  fc  versus  n  . 

Dodge  and  Metzner  (5)  have  performed  a  theoretical 
analysis  for  the  turbulent  flow  of  power-law  fluids  in  smooth  tubes. 

The  following  relationship  was  proposed  between  the  pressure  drop  and 
the  mean  flow  rate, 

(?)  =  Aln  los  (Re'  M1  +  cn'  . (^!) 


where  C  .  =  An  log 

n  In  0 


1  f6n’  +  2  \n< 

"F  V - 5* - ) 


+  c 


n 


the  constants  A^n  and  Cni  were  found  experimentally  and  Equation  4l 
finally  arranged  in  the  form 

vV2  „  /  1-  £*\ 


(?)  =  T^P-75  (f) 


°  4  1 

W)1'2 


(42) 


Equation  42  is  solved  graphically  and  presented  in  the  form 

Dn'  v2_n'  o 


friction  factor  versus  power-law  Reynolds  number. 
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this  plot  theoretical  lines  obtained  for  different  values  of  the  exponent 
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n*  fall  parallel  to  the  Karman-Nikuradse  line  (n1  =  l).  For  values 
of  n*  greater  than  unity  the  lines  are  above  the  n*  =1  line  while 
for  values  of  n*  less  than  unity  they  are  below. 

A  similar  method  of  presentation  has  been  used  by  Shaver 
and  Merrill  (23).  For  experimental,  turbulent  data  a  plot  was  made 
of  friction  factor  versus  power-law  Reynolds  number;  this  yielded 
a  family  of  lines  for  values  of  n'  between  0.53  and  unity.  These 
lines  were  not  parallel,  as  proposed  by  Dodge  and  Metzner  (5),  but 
had  increasing  slopes  for  decreasing  values  of  n',  the  limiting  case 
being  45°  or  an  extension  of  the  Boiseuille  line. 

Ill  Determination  of  Rheological  Properties  Using  a  Rotational  Viscometer. 

A  rotational  viscometer  consists  of  an  inner  and  outer 
cylinder  so  arranged  that  the  fluid  to  be  tested  fills  the  annular 
space  between  the  cylinders.  One  of  the  cylinders  may  be  rotated, 
causing  a  torque  on  the  other  cylinder  which  can  be  measured  and  taken 
as  an  indication  of  the  consistency  of  the  fluid. 

For  a  Newtonian  fluid  the  relationship  between  the  speed  of 
the  rotor  and  the  rate  of  shear  at  the  wall  of  the  stator  may  be 
determined  mathematically  from  the  Reiner-Riwlin  (21)  equation 


where  j  =  rate  of  shear  at  r,  sec"-*-,  when  R;i_<sr<^R2 


(43) 
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N  =  rotor  speed,  revolutions  per  second 
R-^  =  inner  radius  of  the  annulus 
R2  =  outer  radius  of  the  annulus 
If  the  outer  cylinder  is  rotated,  the  rate  of  shear  at 
the  wall  of  the  inner  cylinder  is  required,  so  substituting  Rj_  for  r 

I  dV')  n  R22  N _ 4  n  N 

<^r 'R]_  R2^  -  Rj2  1  -  1/S2 


m 


R2 

where  S  =  — 

R1 

A  linear  relationship  exists  between  the  torque  on  the  inner 
cylinder,  and  the  shear  stress  at  the  wall  of  this  cylinder.  This 
relationship  is 

\  =  K2  9  . (45) 

the  constant  is  obtained  by  calibration  of  the  instrument  with  a 
Newtonian  fluid  of  known  viscosity  at  a  given  temperature . 

Maintenance  of  laminar  flow  in  the  annulus  is  essential.  As 
a  criterion  for  the  onset  of  turbulent  flow  Merrill  (15)>  quoting 
Barr  (2),  suggested 

2  jt  R2  (R2  “  Rj_)  P  N 
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592 


.  .  (46) 


. 

-  -  ■ 

•  ■  •  •'  -  ■'  •  - :  /  .  ... 


n 


•  ■ 


•  •  . 


20 


where 


2  jt  Rg  (R2  ~  Rj)  p  ^ 


is  the  modified  Reynolds  number  of  Barr. 


Krieger  and  Maron  (13)  proposed  a  more  general  equation,  the 
validity  of  which  was  independent  of  the  shear  stress  -  shear  rate 
relationship  for  the  fluid.  Hence  an  expression  for  the  rate  of  shear 
at  the  stator  wall  can  be  obtained.  This  equation  applied  to  the  wall 
of  the  inner  cylinder  is 


where 


it  N 


--  l/s‘ 


1  +  k-L 

S2  - 

2S2 


2 

3 


In 


(W 

(48) 


2 

k  =  s  .  In  S  . (49) 

2  6s2 

n”  =  slope  of  the  logarithmic  plot  of  torque  versus 
rotational  speed. 

Often  the  last  term  of  Equation  47  is  relatively  small  so 
that  the  equation  can  be  reduced  to 


4  jt  N 

1  -  1/s2 


1  +  kx 


.  .  (50) 


Equation  50  can  be  applied  to  either  a  Newtonian  or  a  non-Newtonian  fluid. 

A  paper  published  recently  by  Govier  (9)  discusses  the 
application  of  these  rheological  equations  to  data  obtained  for 
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Bingham  -  and  pseudoplastics  from  a  concentric  cylinder  viscometer. 
The  interpretation  of  the  resulting  shear  stress  -  shear  rate  curves 
for  prediction  of  pressure  drops  in  pipeline  flow  is  discussed  using 
some  of  the  theory  dealt  with  in  section  II. 
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EXPERIMENTAL  EQUIPMENT 

The  Pipeline  Unit 

A  diagram  of  the  flow  circuit  is  shown  in  Figure  3*  The 
flow  circuit  consisted  of  a  supply  tank,  a  pump  and  motor,  the  two 
test  pipelines,  a  venturi  meter,  and  a  weight  tank  on  scales. 

I  Supply  Tank : 

An  18  cubic  foot  capacity  tank  was  used  to  store  the 
suspension,  this  was  connected  directly  to  the  pump  by  a  length  of 
2-l/2  inch  diameter  pipe.  To  ensure  uniform  mixing  of  the  suspension 
two  Greey  "Lightning”  mixers,  one  with  a  paddle  type  impeller  and  the 
other  with  a  propellor  type  impeller,  were  used. 

II  Pump  and  Motor: 

The  pump  used  was  a  "Moyno"  positive  displacement  unit 
operating  at  770  R.P.M.  to  give  35  g.p.m.  against  a  maximum  pressure  of 
75  psi.  The  drive  unit  was  a  5  HP.,  three  phase,  220  volt  motor. 

The  two  basic  pumping  elements  of  a  "Moyno"  pump  are  a  helical 
rotor  and  a  double  thread,  rubber  lined,  helical  stator.  The  two 
helical  surfaces  mesh  together  to  push  the  liquid  ahead  with  uniform 
motion. 

III  Pipeline  Test  Sections: 

The  discharge  from  the  pump  was  designed  to  deliver  suspension 
either  to  the  pipeline  or  return  it  to  the  supply  tank.  Both  this  return 
and  the  discharge  from  the  pipeline  entered  the  supply  tank  tangentially 
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and  well  below  the  surface  of  the  suspension. 

Two  lengths  of  1-7/8  I.D.  aluminum  tube  were  installed 
parallel  to  one  another  and  connected  in  series  by  a  l80°  bend. 

An  iron-constantan  thermocouple  was  inserted  into  the  pipeline  at 
the  bend  to  measure  the  temperature  of  the  flowing  fluid. 

The  pressure  tappings  on  the  section  upstream  of  the  bend 
were  spaced  45*35  feet  apart  and  those  on  the  downstream  section 
spaced  40.17  feet  apart. 

IV  Venturi  Meter: 

The  venturi  meter  was  specially  machined  from  a  cylinder 
of  aluminum  to  give  a  pressure  drop  of  about  15  psi  at  the  maximum 
flow  rate,  and  was  positioned  at  the  end  of  the  downstream  section 
just  before  the  line  entered  the  supply  tank. 

V  Weighing  Tank: 

The  weighing  tank  was  mounted  on  a  set  of  Fairbanks  scales 
with  a  capacity  of  500  lbs.  These  scales  were  fitted  with  a  weigh 
arm  graduated  in  half  pounds.  Emptying  of  the  tank  was  effected  by 
direct  pumping  into  the  supply  tank.  The  line  connecting  the  tank 
with  the  pump  was  initially  a  jointed  metal  hose,  however  this  corroded 
badly  and  was  replaced  by  a  rubber  hose. 

VI  Sample  Cock : 

The  sample  cock  S  (Figure  3)  was  placed  on  the  discharge 
side  of  the  pump  where  a  representative  sample  of  flowing  material 
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could  be  obtained. 

VII  Pressure  Measuring  System: 

A  Statham  electrical  differential  pressure  (DP)  cell 
(Figure  4),  having  a  working  range  of  -  15  psid  was  used  to  measure 
the  pressure  drop  across  each  of  the  test  sections  and  across  the 
venturi  meter.  Figure  5  is  a  schematic  illustration  of  the  electrical 
circuits  associated  with  this  DP  cell.  In  Figure  5  tbe  left  hand  side 
represents  the  excitation  circuit  which  was  made  up  of  an  18  volt  DC 
source,  a  variable  resistance  and  a  voltmeter.  The  right  hand  side 
of  this  Figure  represents  the  output  circuit  which  was  made  up  of  a 
'backed-up'  potentiometer  and  a  recording  potentiometer. 

In  order  to  transmit  the  pressure  to  the  DP  cell  hydraulic 
connecting  lines  were  installed.  These  lines  joined  the  pressure 
tappings  to  one  of  two  four-way  valves  (Figure  6),  depending  on 
whether  it  was  an  upstream  or  downstream  pressure  measurement.  A 
single  line  connected  each  valve  to  the  appropriate  side  of  the 
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Fig.  4.  Differential  Pressure  Cell 
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Fig.  6 


Four-way  Valve  Arrangement 
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EXPERIMENTAL  PROCEDURE 

A.  Initial  Treatment  of  the  Dry  Solids. 

Six  tons  of  coarse  lump  sulphur  were  obtained,  from  Shell 
Oil  Co.  of  Canada  Ltd.,  Calgary.  A  sieve  size  analysis  carried  out 
on  a  representative  sample  of  this  material  showed  it  to  contain 
particles  varying  from  2  inches  in  diameter  to  less  than  37  microns 
(400  mesh).  The  larger  lumps  -were  removed  by  passing  all  the  material 
over  an  8  mesh  “Dillon"  vibrating  screen.  This  large  material  was 
then  crushed,  using  a  heavy  roller,  and  rescreened. 

During  subsequent  screening  an  attempt  was  made  to  obtain 
a  good  separation  of  the  material  into  definite  size  ranges.  The  size 
ranges  sought  had  the  following  Tyler  mesh  values:  8-10,  10-14,  14-20, 
20-28,  28-35 »  35-48,  48-65,  65-100  and  less  than  100.  Due  to  the 
limitations  of  the  screening  equipment  however,  these  size  groups  were 
not  obtained. 

Two  suspensions  were  used  in  this  investigation.  The  first 
was  composed  of  material  obtained  by  the  above  screening,  in  the  size 
range  48-100  mesh,  supplemented  by  larger  material  in  the  range  8-48 
mesh.  Due  to  inefficient  screening  and  lengthy  periods  of  pumping  this 
suspension  (No.  l)  had  the  size  analysis  shown  in  Figure  19  when 
satisfactory  data  were  obtained.  Suspension  "No.  2"  was  made  up  of  a 
certain  weight  of  suspension  1  and  350  lbs.  of  sulphur  flour.  The 
size  analysis  for  this  is  also  shown  in  Figure  19. 
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B.  Operation  of  the  Pipeline. 

The  procedure  adhered  to  from  start-up  of  the  pilot  pipeline 
■was  as  follows: 

1.  The  "Lightning"  mixers  were  switched  on  and  the  suspension 
agitated  until  a  uniform  concentration  had  been  attained. 

2.  The  pump  was  switched  on  and  the  suspension  circulated 
through  the  pipeline. 

3.  Standardization  of  both  the  potentiometers  was  effected 
and  the  D.C.  voltage  applied  to  the  DP  cell. 

4.  Water  from  the  main  source  was  flushed  through  the 
hydraulic  lead  lines  and  the  DP  cell  to  remove  air  bubbles  and  any 
suspension  present  in  the  pressure  tappings. 

5.  The  four-way  valves  were  manipulated  so  that  the  pressure 
drop  over  the  venturi  meter  could  be  measured.  This  pressure  drop 
caused  the  DP  cell  to  produce  a  given  voltage,  which  was  partially 
balanced  by  imposing  a  reverse  E.M.F.  from  the  standard  potentiometer. 
The  residual  voltage  was  then  measured  by  the  recording  potentiometer. 

6.  Manipulation  of  the  four-way  valves  enabled  the  pressure 
drops  over  both  of  the  test  sections  to  be  measured  as  described  under 
section  5» 

7.  The  flow  conditions  in  the  pipeline  were  changed  by  the 
adjustment  of  valves  1  and  2  (Figure  3)  and  the  procedure  for  obtaining 
the  pressure  drops  (sections  5  and  6)  repeated.  Conditions  for  a 


number  of  flow  rates  were  measured  in  this  manner. 
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When  a  complete  series  of  pressure  drop  measurements  had 
been  made  three  samples  were  taken.  These  were 

a)  A  calibrated,  two  gallon  bucket  which  gave  a  bulk 
density  test.  This  sample  was  returned  to  the  supply  tank. 

b)  Four  1  litre  Erlermeyer  flasks  which  were  stoppered 
and  kept  for  subsequent  rheological,  percent  solids  and  size 
distribution  tests. 

c)  A  litre  graduated  cylinder  from  which  settling  data 
was  obtained. 

C.  Testing  Procedure  on  the  Samples  of  Suspension. 

Subsequent  tests  were  carried  out  on  the  suspensions  to 
determine  the  following  properties: 

i)  density 

ii)  initial  settling  rate 

iii)  weight  percent  solids 

iv)  rheology 

v)  particle  size  distribution 

vi)  particle  shape. 

These  tests  are  described  below 

i)  Density: 

The  bucket  containing  sample  (a)  was  weighed  and,  as  the 
volume  was  known,  the  density  could  be  obtained  directly. 

ii)  Initial  Settling  Rate: 

A  stopwatch  was  used  to  determine  the  rate  of  fall  of  the 
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clear  water  suspension  interface  in  a  litre  graduated  cylinder 
( sample  c ) . 

iii)  Weight  Percent  Solids: 

One  of  the  samples  collected  in  a  litre  Erlermeyer  flask  was 
transferred  to  a  beaker  of  known  weight  and  weighed.  The  sample  was 
placed  in  a  drying  oven  and  left  until  no  further  loss  in  weight  was 
observed.  The  weight  of  the  beaker  and  dry  sample  was  then  taken. 

iv)  Rheological  Properties: 

A  concentric  cylinder  rotational  viscometer  was  used  to 
determine  the  rheological  properties  of  the  sulphur-water  suspensions. 

The  particular  model  used,  illustrated  in  Figures  7  and  8,  was  the 

"Fann  V-G  Meter"  (Model  35)  manufactured  by  the  Fann  Instrument  Corporation, 

Houston,  Texas . 

Fluid  is  contained  in  the  annular  space  between  two  coaxial 
cylinders.  The  annular  space  is  filled  to  a  predetermined  height  with 
the  fluid  to  be  tested.  The  outer  cylinder,  or  rotor,  suspended  in  a 
precision  bearing,  is  driven  at  a  constant  rotational  speed  and  the 
torque,  arising  from  the  viscous  drag  of  the  fluid,  is  exerted  on  the 
inner  cylinder,  or  bob.  The  torque  is  rapidly  balanced  by  a  helically 
wound  spring.  The  deflection  of  the  spring,  which  is  proportional  to 
the  torque,  is  read  on  a  calibrated  dial  through  an  optical  reticule. 

In  this  viscometric  method  the  flow  data  are  analysed  in 
terms  of  the  fundamental  relationship  between  shear  rate,  proportional 
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Fig.  7»  Modified  Fann  Viscometer 
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Fig.  8.  Diagram  of  Farm  viscometer 
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to  rotor  speed,  and  shear  stress,  proportional  to  the  dial  reading. 

The  cylindrical  rotor  and  bob  design  permits  calculation  of  the  shear 
rate  corresponding  to  each  rotational  speed.  The  shear  stress  is 
only  obtainable  by  calibration  of  the  instrument  with  a  Newtonian  fluid 
of  known  viscosity.  Careful  calibration  of  the  instrument  was  carried 
out  separately  by  Lindley  (l4)  and  by  Hahn  (ll),  the  calibration  factor 
found  by  these  workers  being  used  in  this  investigation.  This 
calibration  factor  was  checked  using  one  standard  oil.  Hahn  also  found 
that,  at  low  values  of  dial  reading,  a  correction  had  to  be  applied 
to  this  calibration.  The  correction  curve  is  included  in  Appendix  B. 

The  Farm  viscometer  is  a  commercial  instrument;  it  is 
reasonably  accurate  and  simple  to  operate.  Rotor  speed  changes  for  3 > 

6,  100,  200,  300  and  600  R.P.M.  are  effected  smoothly  and  rapidly  without 
stopping  rotation.  This  is  achieved  by  the  use  of  a  three  position  gear 
shift  in  conjunction  with  a  two  position  switch  controlling  the  motor 
speed.  The  rotor  is  driven  by  a  115  volt  50  cycle  dual  speed  synchronous 
motor.  In  the  standard  model  the  sample  being  tested  is  placed  in  a 
stainless  steel  cup,  the  base  of  which  is  equipped  with  pins  to  locate 
and  to  lock  to  a  platform.  This  platform  can  be  moved  vertically  on 
sliders.  The  rotor  and  bob  may  be  detached  from  the  main  body  of  the 
instrument  for  cleaning. 

One  modification  was  made  to  the  standard  model  in  order  to 
allow  further  rheological  data  to  be  obtained.  This  modification  was 
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the  interposing  of  a  9:1  reduction  gear  between  the  synchronous  motor 
and  the  body  of  the  instrument.  This  resulted  in  the  instrument  being 
capable  of  producing  twelve  different  shear  rates  ranging  from 
approximately  1000  reciprocal  seconds  to  l/2  reciprocal  second. 

A  sample  of  the  suspension  was  placed  in  a  600  ml.  beaker 
and  stirred  vigorously  using  an  electric  mixer.  The  beaker  was  placed 
upon  the  movable  platform  of  the  viscometer  and  the  platform  raised  until 
the  rotor  and  enclosed  stator  were  submerged  to  the  depth  indicated  by 
a  mark  etched  on  the  rotor.  The  motor  was  started  and  a  torque,  arising 
from  the  viscous  drag  of  the  fluid  on  the  bob,  was  balanced  by  the 
helically  wound  suspension  spring.  After  the  instrument  was  started  the 
dial  reading  was  observed  to  rise  to  a  maximum.  This  maximum  dial  reading 
and  the  rotor  speed  were  noted.  The  motor  was  stopped,  the  sample  lowered 
and  the  rotor  removed.  Careful  cleaning  of  both  the  rotor  and  the  stator 
was  effected  and  the  instrument  reassembled.  This  procedure  was  repeated 
three  times  for  each  rotor  speed,  an  average  value  of  the  three  deflections 
being  recorded  and  used  to  calculate  the  shear  stress. 

The  number  of  shear  stress  determinations  carried  out  on  a 
given  suspension  was  dependent  upon  the  magnitude  of  the  dial  reading 
at  low  shear  rates.  Consequently  for  suspensions  having  low  percent 
solids  content  the  only  meaningful  data  were  obtained  at  shear  rates  in 
the  range  from  200-1000  reciprocal  seconds,  whereas,  at  high  percent  solids 
concentrations,  readings  could  be  obtained  in  the  range  from  5-1000 
reciprocal  seconds. 
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All  the  rheological  data  given  in  Appendix  B  were 
determined  as  described  above. 

v)  Fhrticle  Size  Distribution: 

A  sample  of  the  suspension  was  weighed  then  wet  sieved, 
using  water,  through  the  following  Tyler  standard  mesh  sieves:  20,  28, 
35j  48,  65,  100,  150,  200,  230,  27O,  325  and  400,  all  of  known  weight. 
The  sieves  and  samples  were  dried  and  then  shaken  on  an  electric  shaker 
to  ensure  a  complete  separation.  Each  sample  was  then  weighed  while 
still  on  the  sieve.  The  weight  of  particles  less  than  400  mesh  was 
obtained  by  difference. 

vi)  Particle  Shape: 

Two  photomicrographs  were  obtained  using  a  Reichart 
Metallograph.  One  was  at  70  magnifications  and  was  of  particles  in 
the  size  range  210-74  microns,  the  other  was  at  140  magnifications  and 
was  of  particles  in  the  size  range  smaller  than  7^  microns. 
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EXPERIMENTAL  RESULTS 


The  experimental  results  are  presented  under  two  headings: 

I  Pipeline  Flow  Data 
II  Properties  of  the  Suspensions 
I  Pipeline  Flow  Data 

The  results  of  the  preliminary  test  on  the  single  phase  flow 
of  water  are  presented  (Figure  9)  in  the  form  of  the  conventional  plot 
of  friction  factor  versus  Reynolds  number.  The  experimental  data 
show  good  agreement  in  the  turbulent  region  with  the  Karman-Nikuradse 
line.  Because  of  the  good  agreement,  it  is  concluded  that  the  cali¬ 
brations  of  the  venturi  meter  and  the  DP  cell  were  accurate,  and 
that  the  auxiliary  equipment  was  functioning  properly.  The  data  for 
this  test  are  tabulated  in  Appendix  C,  Table  1C. 

The  results  of  the  pipeline  flow  tests,  using  sulphur  -  water 
suspensions  1  and  2,  are  presented  in  Figures  10  to  17.  These  are  plots 
of  friction  factor  versus  a  modified  Reynolds  number  (Rew).  This 
modified  Reynolds  number  is  based  upon  the  density  of  the  flowing 
fluid  and  the  viscosity  of  water  at  the  pipeline  conditions.  The 
numerical  data  are  tabulated  in  Appendix  C.  Tables  1C  and  2C. 

Each  of  these  curves  show  two  distinct  portions: 

1.  At  high  Reynolds  numbers  the  curve  is  approximately 
parallel  to  but  displaced  upwards  from  the  Karman-Nikuradse  line  for 
single  phase  fluids.  It  is  apparent  that,  for  both  suspensions,  the 
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separation  "between  the  experimental  and  theoretical  curves  increases 
with  increasing  percent  solids. 

2.  At  lower  Reynolds  numbers  the  data  fall  on  lines  which 
lie  at  an  angle  of  approximately  60°  to  the  abscissa.  For  suspension  1 
these  lines  are  approximately  parallel  and  move  closer  to  the  Poiseuille 
line  as  the  percent  solids  decreases.  For  suspension  2  there  is  a 
definite  tendency  for  the  slopes  to  become  less  as  the  percent  solids 
decreases. 

II  Properties  of  the  Suspensions 

a)  Density: 

The  results  of  the  density  determinations  on  the  suspensions 
are  given  in  Appendix  F,  Table  IF.  These  data  are  plotted  in  Figure  18 
as  specific  gravity  versus  weight  percent  sulphur,  a  comparison  being 
made  with  values  calculated  assuming  a  specific  gravity  of  2.05  for 
sulphur,  zero  solubility  and  additive  volumes.  There  is  satisfactory 
agreement  between  the  experimental  and  calculated  data.  The  small 
displacement  could  be  due  to  air  entrained,  the  presence  of  solid 
impurities  in  the  suspension  or  a  discrepancy  in  the  value  assumed 
for  the  specific  gravity  of  the  sulphur. 

b)  Particle  Size  Distribution: 

Figure  19  is  a  plot,  on  logarithmic  -  probability  co-ordinates, 
of  particle  size  versus  percent  solids  passing  a  given  micron  size. 
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Curves  axe  presented  for  the  two  suspensions  on  which  flow  tests 
were  made  and  for  the  sulphur  flour  used  as  an  ingredient  in  Suspension 
2.  Unfortunately  the  size  analyses  of  the  suspensions  were  terminated 
at  ^00  mesh  and  the  samples  were  destroyed  before  it  was  realized 
that  size  analysis  to  smaller  diameters  was  important.  The  curves, 
representing  the  size  analyses  below  400  mesh  (37  microns)  are 
estimated  but  are  thought  to  be  realistic.  The  size  distribution 
curves  are  of  abnormal  shape  reflecting  the  synthesis  of  the  mixture 
and  the  degradation  by  attrition.  Suspension  2  is  coarser  than 
suspension  1  in  the  first  40  percent  fraction,  slightly  finer  in  the 
next  50  percent  fraction  and  coarser  in  the  last  10  percent  fraction, 
c)  Particle  Shape: 

Figure  20A  is  a  photomicrograph  of  particles  in  the  size 
range  65  to  200  mesh  for  suspension  2,  the  separation  being  obtained 
by  wet  sieving  a  sample  of  suspension.  The  uniformity  in  the  size 
of  the  particles  is  apparent,  this  is  an  indication  that  the  sieving 
procedures  employed  were  efficient. 

Figure  20B  is  the  photomicrograph  of  particles  in  the  size 
range  less  than  200  mesh,  also  for  suspension  2. 

All  particles  appear  to  have  a  fairly  uniform  shape  being 
rounded  with  few  sharp  edges.  This  rounding  of  the  edges  is  probably 
due  to  attrition  during  pumping. 
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Fig.  20 A.  Photomicrograph  of  Sulphur  Particles 
in  Size  Range  65-200  Mesh  (  x  70  Magnifications) 


Fig.  20B.  Photomicrograph  of  Sulphur  Particles 
in  Size  Range  Smaller  Than  200  Mesh  (  x  140  Magnifications) 


~u£o-i  .  is  r  jj.r  ....  '  .  .  v  .  ' 


•  • 


-  53  - 


d)  Initial  Settling  Rates: 

The  initial  settling  rates  of  the  two  suspensions  were 
obtained  by  extrapolating  the  observed  settling  rate  data  to  zero 
time.  Figure  21  shows  this  extrapolated  initial  settling  rate 
plotted  against  weight  percent  solids.  The  data  are  shown  extra¬ 
polated  to  zero  percent  solids,  the  free  falling  velocity  of  the 
largest  particle  is  the  value  used,  this  was  obtained  from  the  Stokes 
equation. 

Settling  rates  for  suspension  1  are  much  less  than  for 
suspension  2  at  the  same  weight  percent  solids.  This  agrees  with 
the  size  analyses  which  show  that  suspension  1  contains  more  fine 
particles  than  suspension  2.  Thus  suspension  1  would  be  expected  to 
have  a  higher  apparent  viscosity.  The  initial  settling  rates  are 
tabulated  in  Appendix  F,  Table  IF. 

e)  Rheological  Properties: 

The  experimental  data  from  the  viscometer  tests  are  tabulated 
in  Appendix  B,  Table  IB.  Initially  these  data  were  plotted  as  dial 
reading  versus  revolutions  per  second.  The  Krieger-Maron  equation 
(Equation  50 )  was  used  to  determine  the  shear  rates,  the  value  of 
the  slope  (n" )  being  found  from  a  smoothed  plot  of  n”  versus  revolutions 
per  second.  The  values  of  the  terms  in  Equation  50  and  the  final  shear 
rates  at  the  wall  of  the  stator  are  tabulated  in  Appendix  A,  Table  1A. 
The  shear  stress  was  obtained  from  Equation  45  using  the  value  of  the 
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calibration  constant  determined,  by  Lindley  (l4)  and  Hahn  (ll),  and 
the  adjustments  to  this  found  by  Hahn. 

Figures  22  and  23  are  logarithmic  plots  of  shear  stress  versus 
shear  rate  for  the  data  tabulated  in  Appendix  B,  Table  IB.  A  small 
constant  has  been  added  to  the  ordinates  and  abscissa  so  that  the  lines 
can  be  extrapolated  to  zero.  Figures  24  and  25  are  arithmetic  plots 
of  the  same  data. 

The  data  indicate  (Figures  22  and  23)  that  both  the  suspensions 
had  complex  rheological  behavior.  Shear  stresses  were  obtained  for 
several  shear  rates  in  the  range  between  70  and  1100  reciprocal  seconds. 
No  indication  of  yield  stresses  were  obtained  during  the  tests  which 
suggested  that  the  materials  were  not  true  Bingham-plastic s.  Figures 
24  and  25,  however,  do  suggest  near  Bingham -plastic  behavior.  Thus 
two  interpretations  are  reasonable,  first  that  of  near  Bingham-plastic 
behavior,  and  second  that  of  complex  pseudoplastic  behavior.  The 
heavy  dashed  lines  on  Figures  22  and  23  indicate  the  pseudoplastic 
interpretation  with  extrapolation  to  lower  and  higher  rates  of  shear. 

Also  shown  on  the  figures  are  light  dashed  lines  which  represent 
"power-law”  pseudoplastic  interpretation  over  a  restricted  range  of 
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INTERPRETATION  OF  THE  PIPELINE  FLOW  DATA 


I  Interpretation  Assuming  Bingham-Plastic  Behavior 

For  interpretation  of  the  flow  data  by  the  nearest  equivalent 
Bingham-plastic  two  physical  properties  were  obtained  from  the  rheo¬ 
logical  data.  These  two  properties  were  the  apparent  coefficient  of 
rigidity  (v^)  and  the  apparent  yield  stress  (Ty&). 

Figure  2 6  is  a  semi-logarithmic  plot  of  apparent  coefficient 
of  rigidity,  in  centipoise,  versus  weight  percent  solids.  The  curves 
are  extrapolated  to  zero  percent  solids  as  the  viscosity  of  water  is 
known  accurately  at  70°F.  The  coefficients  of  rigidity,  at  different 
percent  solids,  were  obtained  by  extrapolation  of  the  rheological  curves 
to  straight  lines  at  shear  rates  greater  than  1200  reciprocal  seconds. 

rn 

In  the  determination  of  xy  the  shear  stress  was  obtained 

a 

for  two  values  of  shear  rate  at  each  percent  solids  concentration. 

Ty0  was  found  by  substitution  of  these  values  of  shear  stress  and  shear 

Cl 

rate  in  Equation  2.  This  resulted  iniairs  of  points,  these  are  shown 
joined  by  dotted  lines  in  Figure  27  which  is  a  semi-logarithmic  plot 
of  apparent  yield  stress  versus  weight  percent  solids.  The  extrapolation 
to  zero  percent  solids  is  again  indicated  as  the  yield  stress  of  water 
is  zero. 

Using  the  values  ofr^a  and  T  ya  taken  from  Figures  2 6  and  27 
the  values  of  the  dimensionless  ratios  Rmg  (Equation  18)  and  Ya 
(Equation  19),  the  apparent  Bingham  Reynolds  number  and  the  apparent 


.^>  --G  .C _ ^  o 


'(.  .j  ’  j 


o_ 


■■  .) 


o  _ _ _ , 


(  ) 


_ 


j 


c  ;  U- 


Q 


_VJ  .  -  O 


VJ 


_  J_ 


L> 


..  LL 


'oJ; 


-j 


.»o 


.  w'  . 


:gu'.I  C 


vy 


-j  : 


■» 


u 


O'  %j  Li 


o 


1\ 


~± 


L 


) ., 


.( 


70 


-  6l  - 


/ 

4-  si 

LEGEND 

ISPENSION 

ISPENSION 

i 

i 

1  O  U 

o  SL 

1 

2 

M 

o 

/ 

’  / 

O 

A 

/ 

V  / 

/ 

X 

1  S  1 

<'+  P 

/ 

y 

o 

/ 

4 

Fig. 

y 

26.  Rheolog 

ical  Data 

Apparent  c 

Defficient  o 

f  rigidity  v 

p.  Weight  pe 

rcent  solids 

20  30  40  50 


WEIGHT  PERCENT  SOLIDS 


o 


10 


60 


70 


A  n  rA  ntiM  i  r  it  lu  iv  u.i ,  aynes  per^rn. 


0  10  20  30  40  50  60  70 


WEIGHT  PERCENT  SOLIDS 


-  63  - 


Bingham  yield  number  respectively,  could  be  calculated  for  experimental 
data  at  different  percent  solids.  Figures  28  and  29  are  logarithmic 
plots  of  the  product  fR^  versus  YQ  for  laminar  flow  data  obtained 

1  Q.  <A 

for  suspension  1  and  suspension  2  respectively.  The  solid  line  in  each 
figure  indicates  the  form  of  the  curve  required  by  Equation  24. 

Figure  28  indicates  that,  by  approximating  the  rheological 
behavior  of  suspension  1  to  that  of  an  equivalent  Bingham-plastic, 
near  agreement  is  obtained  between  the  experimental  and  theoretical 
values  of  fl^  and  Y.  A  similar  approximation  applied  to  suspension 
2  does  not  give  such  near  agreement  (Figure  29).  The  difference  in 
these  results  can  be  attributed  to  the  different  rheological  properties 
of  the  suspensions;  this  difference  has  been  discussed  previously. 

The  data  for  both  suspensions  are  displaced  upwards  from  the  theoretical 
line  indicating  that  neither  are  true  Bingham-plastic s.  It  should  be 
noted  that  the  higher  the  percent  solids  concentration  the  better  is 
the  agreement  between  the  experimental  and  theoretical  data. 

Figures  30  and  31  are  logarithmic  plots  of  friction  factor 
versus  apparent  Bingham  Reynolds  number  with  the  apparent  Bingham 
yield  number  as  parameter.  In  these  figures  both  laminar  and  turbulent 
flow  data  are  plotted.  The  comparison  between  the  laminar  flow  data 
and  Equation  24  is  made  by  imposing  lines  representing  values  of  Y 
equal  to  0,  10,  20,  30  and  40  on  the  plot  showing  the  experimental 
data.  It  is  suggested  that  it  would  be  permissible  to  extend  the  curve 


fV  I  /^\ 


-  6h  - 


APPARENT  BINGHAM  YIELD  NUMBER,  YQ 


o 


FRICTION  FACTOR  -  APPARENT  BINGHAM  REYNOLDS  NUMBER  PRODUCT,  fxR 


O  -  65  - 

6 


APPARENT  BINGHAM  YIELD  NUMBER,  YQ 


APPARENT  BINGHAM  REYNOLDS  NUMBER,  DV/9 


-  68  - 


for  Y  =  10  as  a  loop  to  pass  through  the  points  labelled,  in  Figure  30. 
The  turbulent  flow  data  are  compared  with  the  Karman-Nikuradse  line. 

The  comparison  between  the  experimental  and  theoretical 
data  for  the  laminar  flow  region  again  indicates  that  suspension  1 
behaves  more  nearly  as  a  Bingham-plastic  than  does  suspension  2.  The 
experimental,  turbulent  flow  data  for  both  suspensions  compares  well 
with  the  Karman  -  Nikuradse  line. 

II  Interpretation  Assuming  Power-Law  Pseudoplastic  Behavior 

In  this  interpretation  a  logarithmic  plot  of  friction  factor 
versus  power  law  Reynolds  number  is  used  to  compare  the  experimental 
data  with  theoretical  predictions. 

The  power-law  Reynolds  number,  as  defined  by  Equation  34  is, 
n*  2-n* 


Re* 


D  V  P 


(34) 


where  y 


gcK’  8 


n’-l 


and  K* 


K  ( 


3n»  +  1 

4n* 


(35) 

(32) 


The  values  of  K  and  n*  are  determined  from  the  rheological  data, 
plotted  in  Figures  22  and  23,  in  the  following  manner 

a)  The  maximum  wall  shear  stress  is  determined  from  the 
pressure  drop  data. 

b)  The  range  of  shear  stresses  affecting  the  flow  is 
taken  to  be  between  the  limits  of  this  maximum  shear  stress  and  25$ 
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of  this  shear  stress. 

c)  A  straight  line  is  drawn  through  the  rheological  curve 
"between  these  limits. 

d)  The  slope  of  this  line  is  n1 ,  and  K  is  given  by  the 
intercept  at  a  shear  rate  of  unity. 

The  lines  used  in  the  interpretation  are  shown  lightly 
dashed  in  Figures  22  and  23,  and  the  values  of  n*  and  K  are  plotted 
against  weight  percent  solids  in  Figures  32  and  33 •  The  smoothed 
values  of  n*  and  K  are  used  to  determine  K*  from  Equation  32  and  ^ 
from  Equation  35* 

Figures  3^  and  35  are  logarithmic  plots  of  friction  factor 
versus  the  power-law  Reynolds  number  obtained  from  Equation  3^* 

The  laminar  flow  data  are  compared  with  the  power-law  equivalent  of 
the  Poiseuille  line  ( f  =  l6/R(  t)  proposed  by  Metzner  (l7)*  The 
turbulent  data  is  compared  with  the  lines,  proposed  by  Dodge  and 
Metzner  (5),  representing  the  theoretical  position  of  data  for  n* 
equal  to  1.0,  0.6  and  0.3*  The  poiseuille  line  and  the  lines  re¬ 
presenting  these  values  of  n*  are  shown  on  Figures  3^  and  35* 

The  laminar  flow  data  for  suspension  1  (Figure  3^)  are 
displaced  upwards  from  the  Poiseuille  line.  They  lie  parallel  to 
this  line  and  have  only  a  small  variation  over  the  whole  range  of 
solids  concentrations.  The  turbulent  flow  data  are  also  displaced 
upwards  from  the  lines  predicted  by  Dodge  and  Metzner  (5)  but  the 
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data  again  fall  parallel  to  these  lines  and  lie  in  the  expected 
sequence  (the  highest  friction  factors  associated  with  the  highest 
values  of  n*  at  constant  Re*).  For  suspension  2  (Figure  35)  the 
laminar  flow  data  are  very  scattered  and  are  removed  from  the  Poiseuille 
line.  In  general  the  turbulent  flow  data  fall  close  to  the  positions 
predicted  by  Metzner  and  Reid  for  various  values  of  n* .  The  data 
for  the  54.4  percent  suspension,  however,  is  removed  from  the 
predicted  position.  From  these  data  it  can  be  concluded  that  neither 
suspension  1  nor  suspension  2  behaved  as  a  power-law  pseudoplastic. 

Figure  36  is  a  plot  of  the  critical  friction  factor  versus 
the  power  law  expoent  n* .  The  curve  shown  is  that  predicted  by 
Ryan  and  Johnson  (22)  for  power-law  pseudoplastics. 

The  interpretations  discussed  in  this  section  show  that 
there  are  large  deviations  between  the  experimental  and  predicted 
data.  These  deviations  are  not  necessarily  due  to  the  complexity 
of  the  rheological  behavior.  The  following  analysis  of  the  sources 
of  error  in  the  experimental  results  is  designed  to  indicate  explana¬ 
tions  for  these  deviations. 

As  stated  previously  (page  38)  the  pipeline,  and  associated 
pressure  measuring  equipment,  were  calibrated  by  pumping  water  through 
the  equipment.  Good  agreement  was  obtained  between  the  measured  data, 
plotted  as  friction  factor  versus  Reynolds  number  (Fig.  9)  and  the 
Karman-Nikuradse  line.  This  indicated  that  the  equipment  was  function¬ 
ing  satisfactorily  and  that  the  calibration  data  for  the  pressure 
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measuring  system  were  accurate.  The  pressure  drop  data  are,  therefore, 
assumed  to  he  accurate.  This  is  substantiated  by  the  good  agreement 
between  the  experimental  data,  when  plotted  as  friction  factor  versus 
modified  Reynolds  number  (Fig.  10  through  IT)  and  the  Karman-Nikuradse 
line. 

The  possibility  of  errors  in  the  rheological  data,  however, 
is  much  greater.  A  majority  of  these  data  were  obtained  in  the  low 
dial  reading  range  of  the  Farm  V-G  meter  where  the  calibration  curve 
is  no  longer  linear  and  small  errors  in  the  dial  reading  are  magni¬ 
fied.  Thus  errors  in  Ty,  and  n  may  be  large.  No  method  of  deter¬ 
mining  the  magnitude  of  these  errors  is  available. 
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CONCLUSIONS 

In  the  conclusions  the  two  suspensions  will  he  considered 
separately  as  their  rheological  behaviour  is  quite  different. 

1.  Suspension  1 

a)  The  rheological  behaviour  of  this  suspension,  in  the 
range  of  the  pipeline  tests,  resembles  that  for  a  Bingham-plastic.  Thus 

i)  Data  obtained  in  the  laminar  flow  region 
correlate  approximately  with  that  predicted  for  Bingham-plastic s. 

ii)  Friction  factors  for  the  turbulent  flow  region 
are  independent  of  the  apparent  Bingham  yield  number. 

iii)  Data  obtained  in  the  transition  between  laminar 
and  turbulent  flow  follow  a  fairly  well  defined  curve  at  constant 
values  of  the  apparent  Bingham  yield  number. 

b)  The  rheological  behaviour  also  resembles  that  for  a 
power-law  pseudoplastic.  The  agreement  between  experimental  and 
predicted  data  is  not  however  as  close  as  that  obtained  when  Bingham- 
plastic  behaviour  is  assumed.  Thus 

i)  The  laminar  flow  data  correlate  approximately 
with  the  curve  f  =  16/Re 1 . 

ii)  The  turbulent  flow  data  do  not  follow  the 
theoretical  predictions. 

c)  This  suspension  behaves  as  a  pseudoplastic,  however, 
in  the  range  of  the  experimental  data,  it  can  be  treated  as  a 
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Bingham-plastic. 

2.  Suspension  2 

a)  The  rheological  behaviour  of  this  suspension  indicates 
that  it  is  not  a  Bingham-plastic.  Thus  it  is  concluded  that 

i)  The  laminar  flow  data  does  not  correlate  as  a 
Bingham-plastic  at  constant  values  of  apparent  Bingham  yield  number. 

ii)  Data  obtained  in  the  turbulent  flow  region  still 
correlate  approximately  with  the  Karman-Nikuradse  line  as  it  is  independent 
of  the  apparent  Bingham  yield  number. 

b)  As  indicated  by  the  rheological  data  this  suspension 
behaves  as  a  pseudoplastic.  However  it  is  not  a  power-law  pseudoplastic 
as 

i)  The  laminar  flow  data  do  not  correlate  with  the 

curve  f =l6/Re 1 . 

ii)  The  turbulent  flow  data  only  correlate  approximately 
with  the  theoretical  predictions. 
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Nomenclature 

A1 

Constant  defined  by  Equation  4l,  dimensionless. 

c  ,  c  , 

Constants  defined  by  Equation  4l,  dimensionless. 

dV/dr 

Shear  rate,  sec."1  ( dV/dr )_  refers  to  the  shear  rate  at 

til 

the  vail  of  the  inner  cylinder  of  a  rotational  viscometer. 

D 

Diameter  of  tube,  ft. 

E&P 

4l 

Shear  stress  at  vail  of  a  tube,  lb^/ft.2 

f 

Fanning  friction  factor,  dimensionless. 

F 

Shearing  force,  lbp. 

Sc 

p 

Dimensional  conversion  factor,  32.2  ft.  Ib^  sec.  lbp. 

H 

_2 

Ty  D  p  gQ 

Hedstrom  number  x  ,  dimensionless. 

V  k2 

Instrument  constants  defined  by  Equations  48  and  49, 

dimensionless. 

K 

Pover-lav  constant  defined  by  Equation  3  j  lbp/sec.n  /ft.2 

L 

Length  of  tube,  ft. 

n 

Exponent  defined  by  Equation  3  >  dimensionless. 

1 

n 

Exponent  defined  by  Equation  30 »  dimensionless. 

IT 

n 

■  ^  ln  ^  ,  dimensionless, 
d  In  N 

N 

Rotational  speed,  rev. /sec. 

P 

Pressure,  lb^/ft.^  dP^refers  to  pressure  loss  due  to  friction 

Q 

Volumetric  rate  of  flov,  ft3/sec. 
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Re 

Re' 

*m 

RiR2 


r 


s 

T 


t 


Umax 

V 

V* 

8v/d 


v 

W 

Y 

Z 


Newtonian  Reynolds  number,  dimensionless. 

Power-law  Reynolds  number,  dimensionless. 

Bingham  Reynolds  number,  dimensionless. 

Radius  of  the  inner  and  outer  cylinder  in  a  concentric 
cylinder  viscometer,  ft. 

Radial  distance  from  tube  axis,  ft.  rv  is  the  radial 
distance  to  the  tube  wall . 

Ratio  R-[_/R2>  dimensionless. 

Shear  stress,  dynes/cm^.  Tw  refers  to  the  shear  stress 
at  the  wall  of  a  round  tube,  Tp^  to  the  shear  stress  at 
the  wall  of  the  innder  cylinder  of  a  rotational  viscometer, 
and  Ty  is  the  yield  stress  of  a  Bingham-plastic. 

Torque  on  the  spring  of  a  rotational  viscometer. 

Maximum  velocity  of  flow,  ft/sec. 

Local  axial  velocity,  ft/sec. 

Mean  velocity  of  flow,  ft. /sec. 

Shear  stress  velocity,  ft/sec. 

Shear  rate  at  the  wall  of  a  tube  for  a  Newtonian  fluid  in 

laminar  flow,  sec.  1 

Specific  volume,  ft^/lb  . 

Flow  function  Q/  r3  sec."1 

■w 

Yield  number  D  Ty  gc/V  y|^  ,  dimensionless. 

Criterion  defined  by  Equation  36  .  Z^ax  refers  to  maximum 
value  of  this  criterion. 


-  81  - 


z 


y 


9 

©  ’  ,  9",  9"’ 

M- 

P 

0  (n) 

T 


Elevation  above  a  given  datum,  ft. 

Constant  in  Equation  4,  dimensionless. 

Denominator  in  power-lav  Reynolds  number  (Equation  3^)> 
lbm/ft  sec.2“n 
Pipe  roughness,  ft. 

Coefficient  of  rigidity  of  a  Bingham-plastic  fluid,  dynes 
sec./cm^. 

Measurement  of  shear  stress  in  a  rotational  viscometer,  sec."-1- 
Indicators  of  functional  relationships. 

Viscosity  of  a  Newtonian  fluid,  dynes  sec./'cm^  (poise). 

Fluid  density,  lbm/ft3. 

Constant  defined  by  Equation  39 »  dimensionless. 

Shear  stress,  dynes/cm  .  Ty  refers  to  the  shear  stress 
at  the  wall  of  a  round  tube,  to  the  shear  stress  at 
the  wall  of  the  inner  cylinder  of  a  rotational  viscometer, 
and  t  is  the  yield  stress  of  a  Bingham-plastic. 
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APPENDIX  A 

Computation  of  the  shear  rate  in  the  Viscometer 

Supplementary  Data 
Stator  radius  R^  =  O.679  inches 
Rotor  radius  R2  =  0.724  inches 
^  =  5  =  I.O663 

The  shear  rate  at  the  stator  wall  for  a  non-Newtonian  fluid,  is: 
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APPENDIX  B 

Viscometer  experimental  data 

The  shear  stress  at  the  stator  wall  is 

p 

T.,.,  =  Kn  9  dynes  per  cm 

R1  1 

or  for  a  Newtonian  fluid 

Tt,  =  p  (—)  dynes  per  cm^ 
^1  'dr/j^ 

where  the  viscosity  p  is  in  poise. 


It  follows  that 

v  ■  **  m 


However  for  a  Newtonian  fluid 

dV' 

>dr  / 


(£1)  =  K0  N  =  104.3  N 

Vdr/R.  * 


so 

or 


KjQ  =  p  Kg  N 
-  M-  k2  N 


K-, 


9 


The  value  of  determined  by  Lindley ^ ^ ^AND  Hahn was  4.96 


The  correction  curve  for  this  calibration  constant  is  shewn  in  Figure  13. 
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APPENDIX  C 

PIPELINE  EXPERIMENTAL  DATA 
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The  data  tabulated  in  Tables  2G  and  3G  are  shown  plotted  in  Figure  2G  as  flow  rate 
(lb. /hr.)  versus  / AP.P . 


n 


c  • 

r. 


* 


-  _ 


i 

f.  .  -i 


-  .  r  j 


.  T 


■  -'i 


l  .•  '  ^ 

-■  1 


I*  '.;  J 

"  .  -  , 

•  i 


,  .  •> 


*  •  «  • 


I  I 


. 

i  ; 

:•  -i  t 

*  •  • 

•  « 

• 

i ' 

-  \  • 

'  ■) 

• 

• 

